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Abstract 

We present a review of cosmological nucleosynthesis (CN) with neutrino oscil- 
lations, discussing the different effects of oscillations on CN, namely: increase of 
the effective degrees of freedom during CN, spectrum distortion of the oscillating 
neutrinos, neutrino number density depletion, and growth of neutrino-antineutrino 
asymmetry due to active-sterile oscillations. We discuss the importance of these 
effects for the primordial yield of helium-4. 

Primordially produced 4 He value is obtained in a selfconsistent study of the 
nucleons and the oscillating neutrinos. The effects of spectrum distortion, depletion 
and neutrino-antineutrino asymmetry growth on helium-4 production are explicitly 
calculated. 

An update of the cosmological constraints on active-sterile neutrino oscillations 
parameters is presented, giving the values: <5m 2 (sin 2 2$) 4 < 1.5 x 10 -9 eV 2 for 
5m 2 > 0, and |<5m 2 | < 8.2 xlO -10 eV 2 at large mixing angles for dm 2 < 0. According 
to these constraints, besides the active-sterile LMA solution, also the active-sterile 
LOW solution to the solar neutrino problem is almost totally excluded. 



1 Introduction 



Cosmological nucleosynthesis is traditionally used as a probe of the conditions of the 
early Universe at the nucleosynthesis epoch (T ~ MeV). Concerning neutrino physics, 
the requirement for a concordance between the theoretically predicted and the extracted 
from observations primordial abundances of light elements, constrains neutrino charac- 
teristics: mass, number of light species, degeneracy, decay width and mass of eventual 
heavy decaying neutrinos, neutrino oscillation parameters, possible new interactions, etc. 

On the other hand, the topic of neutrino oscillations is with us since Pontecorvo's 
hypothesis for these [|1] i.e. more than 40 years. Neutrino oscillations imply non-zero 
neutrino masses and mixings, therefore presenting an indication of physics beyond the 
electroweak standard model. Recently the positive indications of oscillations obtained 
from the greatest neutrino experiments (SuperKamiokande, SNO, Soudan 2, LSND, etc.) [] 
turned the subject of neutrino oscillations into one of the hottest points of astrophysics 
and neutrino physics. Hence, it looks like appropriate to provide an updated review of 
the influence of neutrino oscillations on CN and present the most recent cosmological 
constraints on neutrino oscillations parameters. 

In case neutrino oscillations are present in the primordial plasma of the early Universe, 
they may lead to changes in CN, depending on the oscillation channels and the way they 
proceed. Namely, the oscillations effect depends on the kind of oscillations (they can be 
resonant or nonresonant) and also differs for the equilibrium and nonequilibrium cases. 
Oscillations may influence nucleosynthesis through their effects on: neutrino and antineu- 
trino number densities, spectrum, neutrino-antineutrino asymmetry and the number of 
neutrino species. 

In the next section we will provide a review of CN with neutrino oscillations, discussing 
mainly oscillations influence on CN. In the last section we will present an update of the 
cosmological constraints on v e v s neutrino oscillation parameters and discuss how they 
concern the solutions to the solar neutrino problem. 



1 All recent data of the solar neutrino experiments namely Gallium, Chlorine, SuperKamiokande, SNO, 
sensitive to different neutrino energy, point to a well pronounced deficit of solar neutrinos in contrast 
with the theoretical predictions of the Standard Solar Model fluxes. 

The measured ratio of the muon over electron neutrino flux by Kamiokande, SuperKamiokande, 1MB 
and Soudan 2 is considerably lower than the predicted one for the atmospheric neutrinos, resultant from 
the primary cosmic rays collisions with the nuclei in the upper atmosphere. Moreover, a zenith-angle- 
dependent deficit of muon neutrinos was observed. 

The LSND experiment has observed electron antineutrino appearance in a flux of muon antineutrinos 
and electron neutrinos in a flux of muon neutrinos. 

These three neutrino anomalies are preferably explained in terms of neutrino oscillations, which are 
able to provide a natural mechanism for energy dependent and neutrino type dependent suppression of 
neutrino fluxes. There exist several solutions to the solar neutrino problem: Small Mixing Angle (SMA), 
Large Mixing Angle (LMA), LOW and Vacuum Oscillations (VO) solutions, which imply neutrino mass 
differences in the range: 10~ 10 — 10~ 4 eV 2 . The atmospheric neutrino anomaly can be resolved by 
oscillations with mass differences 10~ 3 — 10~ 2 cV 2 . While for the LSND neutrino experiment, much 
bigger mass difference, of the order of eV is necessary. 
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2 Historical overview of CN with neutrino oscilla- 
tions 



2.1 Standard Big Bang Nucleosynthesis 

For a precise analysis of the oscillations effect on CN, the element that is used traditionally 
is helium-4, as far as the most reliable and abundant data now available are for that 
element. According to the standard Big Bang Nucleosynthesis (SBBN) theory 4 He is a 
result of a complex network of nuclear reactions, which proceed after the freezing of the 
neutron-to-proton ratio n/p. 

The abundance of the primordially produced mass fraction of helium-4 Y p ~ 2(n/p)f/ 
(1 + (n/p)f) depends mostly on two compelling processes, determining the neutron-to- 
proton freezing ratio - (n/p)/, namely the Universe's cooling rate, H(t) ~ y/9eff T 2 
and the interaction rates of the weak processes, interchanging neutrons and protons: 
T w ~ G 2 F (gy + 3g\) E 2 T 3 . Hence, the produced helium is a strong function of the effective 
number of relativistic degrees of freedom at the CN epoch, g e ff, and the neutron mean 
lifetime r n , which parametrizes the weak interactions strength. Besides, primordially 
produced 4 He is a logarithmic function of the baryon-to-photon ratio r), due to the nuclear 
reactions dependence on nucleon densities, i.e. Y p (g e ff,T n ,T]). Deuterium measurements 
in pristine environments towards low metallicity quasar absorption systems at very high 
z~3 provide us with the most precision determination of the baryon density, giving 
the value: rj = 5.6 ± 0.5 x 10~ 10 . Recently, the baryon density was also determined from 
observations of the anisotropy of the cosmic microwave background (CMB) by DASI, 
BOOMERANG and MAXIMA experiments @. The CMB value is in agreement with the 
one found from deuterium measurements and SBBN. 

SBBN assumes three neutrino flavours, zero lepton asymmetry and equilibrium neu- 
trino number densities and spectrum: 

n e J> = exp(-£/T)/[l + exp(-E/T)] 

The primordial 4 He abundance Y p , predicted from SBBN, is calculated with great preci- 
sion. In recent years, corrections to SBBN code accounting for different physical effects 
reduced the theoretical uncertainty to less than 0.1% (\SY P \ < 0.0002) within a wide range 
of T] H]. The uncertainty of the observational Y p is few percent. The predicted helium-4 
value is in accordance with the contemporary helium values, inferred from observational 
data: 0.238-0.245 (the systematic errors are supposed to be around 0.007) J5], and is 
consistent with other light elements abundances. 

Given this accuracy, it can be used as a probe of the eventual new neutrino physics - 
the neutrino oscillations. 



2 Where T-^G 2 M + 3g 2 A ). 
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2.2 CN with neutrino oscillations 



CN with neutrino oscillations was studied in numerous publications P~P5|. The basic 
idea of oscillations is that neutrinos are not degenerate in mass and there is a mixing, so 
that the mass eigen-states V{ are distinct from the flavour eigen-states Vf and are unitary- 
combinations of the latter: 



Ui = U if u f (/ = e,/i,r). 
Then in the simple two-neutrino oscillation case 



Z/i = CVf + SVf/, 
V 2 = —SPf + CPfi, 

where c = cos("$), s = sin(i9), $ ^ is the mixing angle, v\ and z/ 2 are the mass eigen-states 
with masses correspondingly m\ and m 2 and Sm 2 ^ 0. 

The probability to find at a distance / a given neutrino type /' in an initially homo- 
geneous neutrino beam of another neutrino type / is: 

Pff» = sin 2 2$sin 2 (1.267 5m 2 l/E), 

where Sm 2 is the neutrino mass difference in eV 2 , E is the neutrino energy in MeV and / 
is the distance in meters. 



2.2.1 Vacuum oscillations 

The oscillations effect on CN was historically first considered for vacuum oscillations 
between different neutrino flavours || [7| . It was shown that since there is a slight deviation 
from equilibrium in that case (temperatures of different flavour neutrinos are nearly the 
same), flavour oscillations have a negligible effect on neutrino number densities and on 
nucleosynthesis [|]. The possibility of generation of neutrino-antineutrino asymmetry in 
the electron sector due to CP-violating oscillations and its effect on CN was discussed in 
ref. 0. 

Active-sterile neutrino f] oscillations effect on CN may be noticeably stronger because 
they can 

(a) increase the effective number of light degrees of freedom during nucleosynthesis 0, 



(b) lead to a strong distortion of the neutrino spectrum |10| 



(c) lead to a considerable decrease of the total number density of active neutrinos [|1C 
Active-sterile oscillations may keep sterile neutrinos in thermal equilibrium or bring 
them into equilibrium in case they have already decoupled. CN allows not more than one 
additional neutrino type, therefore forbidding efficient production of sterile neutrinos due 
to oscillations (T r < H). First constraints || on oscillation parameters were obtained 



i SU2 singlet neutrino. 
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from the requirement that the observed primordial values allow less than one additional 
light neutrino type. It was shown that the simultaneous presence of Dirac and Majorana 
neutrino masses would contradict the observed helium-4 abundance if Sm 2 > 10~ 6 eV 2 . 

In ref. || CN constraints were estimated for fast oscillations (T osc > H) from the 
condition that sterile neutrinos should not be abundantly produced due to oscillations 
before nucleosynthesis epoch (Tr ~ T w sm 2 2^ < H at 1 MeV): (Sm 2 ) 2 sin 2 2$ < 1.6 x 
l(T 17 eV 4 . 

In ref. || a condition for an efficient production of sterile neutrinos at about 1 MeV, ac- 
counting for the rescattering of neutrinos from the fermions of the plasma, was calculated 
and constraints on neutrino oscillation parameters were also obtained. These constraints 
H U as will be discussed in more detail in the next subsection, are considerably relaxed, 
when a proper account of the dispersion effects is provided. 

In ref. the nonequilibrium case of active-sterile oscillations, effective after the neu- 
trino decoupling, i.e. the oscillation rate exceeding the expansion rate r osc ~ Sm 2 / (4i?) > 
H, while neutrino weak rates less than expansion rate r„, < H , was discussed. 

The effect of such active-sterile neutrino oscillations is due mainly to the fact that 
sterile neutrinos may have decoupled much earlier than the active ones [] and at the 
nucleosynthesis epoch their number densities were negligible in comparison with the active 
neutrinos: iV s <C Nf, as far as T s < Tf and N ~ T 3 . 

The accurate kinetic approach of ref. H to a description of the oscillating neutrinos 



in terms of neutrino density matrix in momentum space, was used in ref. |lOj to calculate 
analytically the evolution of the number density of electron neutrinos: 

n Ve = p LL = {l - 2cV + 2cV cos[BT/E(T- 3 - T 3 )}} x <«, 

where B = 0.1M P iSm 2 / \fg^f f and T = 3 MeV. 

It was explicitly shown, that for large interval of oscillation parameters of the discussed 
model, oscillations may cause considerable spectrum distortion and/or depletion of the 
electron neutrino. 

At great Sm 2 , cos[BT/ E(T~ 3 — T ~ 3 )] is frequently oscillating and can be averaged. 
Then pll = (1 — 2c 2 s 2 )n^, and active neutrinos number densities are depleted by oscil- 
lations. However, for smaller Sm 2 , the energy distribution of the active neutrinos may be 
also considerably changed. 

Both the electron neutrino depletion and the decrease of the electron neutrino energy 
due to oscillations into less energetic sterile neutrinos decrease the weak rates I 1 ^ ~ 
E 2 N U , leading to higher freezing temperature of the nucleons. The total effect is an 
overproduction of He-4. 

Helium production was numerically calculated for the full range of parameters of this 
oscillation model. It was shown that, the effect of depletion and spectrum distortion 
due to oscillations on He-4 may be much stronger than the effect due to excitation of an 
additional degree of freedom. Assuming that allowed overproduction should be less than 

4 Therefore, they were not heated by the following after their decoupling processes leading to entropy 
increase. 



5 



5%, the following constraints on the oscillation parameters were obtained: 5m 2 < 10~ 9 
at i? > 7r/15 The overproduction of helium due to sterile neutrino species excited by 
oscillations was calculated. 



2.2.2 Matter oscillations 

In the following years, it was realized that for a large interval of masses and mixing 
angles the thermal background in the pre-nucleosynthesis epoch may strongly affect the 
propagation of neutrinos and, has to be taken into account. Differences in the interactions 
with the particles from the plasma lead to different average potentials for different neutrino 
types [|ll|]. In the adiabatic case the effect of the medium can be formally hidden in the 
oscillation parameters, by introducing matter oscillation parameters that are expressed 
through the vacuum ones and through the characteristics of the medium. The matter 



mixing angle is then |11], |39 



sin 2 t? m = sin 2 ^/[sin^ + (g^L-cos2^) i ], 

where Q = -bE 2 T A / '(5m 2 M^) , L = -aET 3 L a /(5m 2 ), L a is expressed through the 
fermion asymmetries of the plasma, a and b are positive constants different for the different 
neutrino types, — L corresponds to the neutrino and +L to the antineutrino case. 

It was realized that although most often the medium suppresses oscillations (by de- 
creasing their amplitude), there also exists a possibility of enhanced oscillation transfer 
in case a resonant condition between the parameters of the medium and the oscillation 
parameters holds: 

Q T L = cos 2tf . 

Then the mixing in matter becomes maximal, independently of the value of the vacuum 
mixing angle. Besides, as can be easily judged from the resonant condition, at high 
temperatures, when \Q\ > \L\, 5m 2 > corresponds to a nonresonant case, while 5m 2 < 
corresponds to a resonant case, and the resonance holds in both neutrino and antineutrino 
sectors. At low temperatures, when \Q\ < |L|, the resonance is possible either for neutrinos 
in the case 5m 2 > or for antineutrinos in the case Sm 2 < 0. 

Matter oscillations of ordinary flavour neutrinos, occurring before the freeze-out of 
n/p - ratio were considered first fl5fl . Two different effects of oscillations were taken into 



account: the generation of v e — v e asymmetry in v T <-> v e oscillations resonant conver- 
sion and the change in the neutrino freeze-out temperature. The effect on primordial 
production of He-4 was estimated to be very low: 5Y P < 1.3 x 1CT 3 . 

The effect of matter active-sterile oscillations on CN, taking into account the essen- 
tial processes of neutrino forward scattering of the background particles, was considered 
first in refs. flIS|,|lI?fl-f2"0fl,[p4l}]. Active-sterile oscillations that are efficient before the nu- 
cleosynthesis epoch (corresponding to high mass differences) were discussed there. The 
production of sterile neutrino states due to oscillations and interactions of the active 
neutrino with the medium was calculated. 

It was shown that matter oscillations, proceeding before the active neutrinos decou- 
pling, are capable of exciting additional degrees of freedom into the plasma. The increase 
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in the expansion rate H ~ y/g e ff and the corresponding overproduction of helium were 
calculated. 

The exclusion regions for the neutrino-mixing parameters were obtained from the 
requirement that the number of neutrino species is less than a certain value (usually N v < 
3.4 in the pioneer publications), up to which there is an agreement between the values 
for light element abundances obtained from theoretical predictions and those extracted 
from observations. The excluded region of oscillation parameters values was found much 
smaller than in the case not accounting for the neutrino coherent interactions. 

In these works the kinetic evolution of the neutrino ensembles was studied in terms 
of particle densities, and it was assumed that oscillations slightly shift neutrino density 
matrix from its diagonal form. First estimates of an eventual depletion of electron neu- 
trinos, due to u e <-> u s oscillations, and its corresponding influence on helium-4 were 
made in [1171 EOfl. Numerical calculations of the discussed oscillations effects on helium-4 



production were made in p4 



The idea that active-sterile oscillations will produce considerable v e — v e asymmetry 



during resonant transfer, which may influence CN was discussed in ref. [16 



2.2.3 CN constraints on active-sterile oscillations 

The pioneer works discussing u e <-> v s at 5m 2 > excluded the Large Mixing Angle 
solution (LMA) of the solar neutrino problem. 

In the Up <->■ u s oscillations case, the solution to the atmospheric anomaly was found 
to lie in the cosmologically excluded region. This can be regarded as an indication that 
the atmospheric neutrino anomaly should be solved by oscillations in the v T sector 

(which is also the experimentally preferred recent solution) []. 

Recently |36| the constraints on active-sterile neutrino oscillations were discussed in 



connection with the controversy in literature concerning the sterile neutrino production 
rate: in some works annihilation rate was used when calculating it, while in others the 
total reaction rate was employed. The problem is still under discussion. 



In ref. [25] (see also ref. [Bq]) a rough estimation of the spectrum distortion by shifting 



the effective temperature of the neutrino and considering its spectrum equilibrium was 



provided. However, as was shown in refs. [29, 30, 31, 33 1 this way the real distortion of 



the neutrino spectrum cannot be described. Besides, the exclusion plots of refs. [|^, £6] 
for the resonant case are not compatible with the ones for the nonresonant case, they do 
not coincide at maximal mixing, as they should. Hence, it is preferable to use the results 
of ref. |23| for the cosmological constraints in the <-> v s and v T v s as more reliable. 

It should be kept in mind also that the constraints on z/^ ]T «-> v s have not been updated: 
they do not take into account the spectrum distortion of the oscillating neutrinos and the 
growth of the neutrino-antineutrino asymmetry. 



5 With this in mind, and also the fact that the asymmetry growth due to resonant oscillations still 
needs some clarification and precision calculations, we will not discuss the numerous works aiming to 
save the active-sterile solution, by weakening the cosmological constraints due to the large asymmetry 
production in active-sterile oscillations in the muon or the tau-neutrino sector. 
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In the resonant case we expect the bounds to become less restrictive in the small mixing 
angles area, when the estimation of the asymmetry growth due to oscillations is made, 
because oscillations generated asymmetry at small mixing angles leads to a suppression 
of oscillations and consequently to a relaxation of the bound on oscillation parameters 
at small mixings []34|| . In both resonant and nonresonant cases we expect less restrictive 
bounds at large mixings due to spectrum distortion effect. 



2.3 Update of oscillations effects on CN 

In the last years the oscillations effects on CN were updated [^J-[[£|,[|37|, by providing 
a precise numerical account for the neutrino spectrum distortion, depletion and growth 
of neutrino-antineutrino asymmetry due to neutrino oscillations. 



2.3.1 Spectral distortion. 



In refs. ||29||-||33||, ||37|| exact kinetic equations for the neutrino density matrix in momen- 
tum space, accounting simultaneously for Universe expansion, neutrino oscillations and 
neutrino forward scattering, were numerically solved for oscillations with small mass dif- 
ferences. The equations for the neutrino density matrix were solved self consistently with 
the kinetic equations of the nucleons during the n/p - freezing (see detail description in 
the next section), which allowed to provide a precise analysis of neutrino depletion and 
spectrum distortion effect on CN due to oscillations. The spectrum distortion effect was 
shown to be considerable fl29j. In Fig. 1 the net effect of spectrum distortion is given by 



the difference between the short-dashed curve from ref. [p4j] and the long-dashed curve 
from ref. [Q. In the nonequilibrium case of electron-sterile oscillations, it cannot be 
estimated simply by shifting the effective temperature or introducing chemical potentials. 

On the other hand the precise numerical analysis of the spectrum distortion is tech- 
nically very difficult. It has been provided only for oscillations with small mass differ- 



ences 



m 



For the adequate description of the spectrum distortion in that case 
1000 bins in the nonresonant case and up to 5000 bins for the resonant case were needed. 
Hence, thousands times larger system of coupled integro-differential equations had to be 
accurately solved, in contrast to the case of momentum averaged calculations when a 
set of 8 ordinary differential equations are solved. The computational time was strongly 
increased. 

Recently a successful analytical description of the spectrum distortion of neutrino 
spectrum due to oscillations was proposed in ref. [|3 for small mixing angles and large 
mass differences. It may be very useful in precision calculations of the oscillations effect 
on CN at large mass differences, which have not been done. 



2.3.2 Neutrino-antineutrino asymmetry growth 

The idea of generating neutrino-antineutrino asymmetry due to resonant transfer of neu- 
trinos in active-sterile neutrino oscillations was first proposed in ref. 
developed further by the authors and remained unnoticed. 



16fl . Alas, it was not 
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Figure 1: Comparison between the results on primordial helium-4 production obtained in our 
works [32, and in previous work The long-dashed curve shows our results when the 



asymmetry effect is neglected. The short-dashed curve shows the results of Enqvist et al., where 
both asymmetry and spectrum distortion effects were neglected. The difference between the two 
curves shows the effect of the proper account of the spectrum distortion of the neutrino. On the 
other hand, the difference between our curves, the solid and the dashed one, presents the net 
asymmetry effect. 

Estimations of the effect of asymmetry, generated in flavor oscillations, on CN of 
helium-4 were provided in ref. p5| |. The effect was shown to be negligible. 

Oscillations generated asymmetry was discussed in more detail for active-sterile os- 
cillations with large mass differences in refs. [17t|-[p0|, where it was estimated to have 



negligible effect on CN, and in ref. |27j, where asymmetry growth was found possible in 
case 6m 2 > 10~ 5 eV 2 . In the numerical study, exploring oscillations with small mass 
differences 5m 2 < 10~ 7 eV 2 | 29| , a growth of asymmetry was registered. 

The revealment of the asymmetry growth at small mass differences became possible 
thanks to the exact kinetic approach, namely work with selfconsistent kinetic equations 
for neutrinos in momentum space. This approach allows to take into account the neutrino 
spectrum distortion and describe the asymmetry behaviour at each momentum. 



In ref. [32] the evolution of the asymmetry was studied selfconsistently with the neu- 
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trino and nucleoli evolution, which allowed to take into account the effect of the asymmetry 
on CN during the full evolution of the asymmetry. It was shown that besides the well- 
known kinetic effect of the asymmetry on CN, after it has grown up to high values of the 
order of ±0.01, very small asymmetries L < 0.01 exert a considerable indirect effect on 
CN through oscillations. This is due to the fact that even very small asymmetries change 
the medium-induced neutrino potential energy and may have a considerable back reaction 



on the oscillating neutrinos. It has been shown |32|, p7| , that the results obtained when 
the indirect asymmetry effect was taken into account differ by many orders of magnitude 
from results neglecting small asymmetries. 

For the case of small mass differences exact kinetic calculations of the asymmetry effect 



were provided in refs. p9|-||34j| . The asymmetry growth up to 5 orders of magnitude was 
registered in these studies. Hence, these lepton asymmetries experienced indirect effect 
on CN through changing neutrino number densities and the oscillation pattern, which on 
its turn effected the kinetics of nucleons during n/p - freezing. 

The net indirect asymmetry effect on CN is illustrated in Fig. 1 by the difference 
between the solid and the long-dashed curves. The asymmetry growth reduced neutrino 
oscillations effect and lead up to 10% underproduction of helium in respect to the model of 
CN with oscillations but neglecting asymmetry. Hence dynamically produced asymmetry 



leads to a noticeable weakening of the cosmological constraints for small mixing angles [29 



However, the precise numerical description of the asymmetry growth requires more 
than 5 000 bins for the spectrum distortion, which was found essential for the correct 



calculation of the asymmetry evolution [ 33j] . This increases enormously the calculational 
time in the resonant case. Besides, the neutrino evolution equations at resonance have 
high stiffness and implicit methods should be used to solve the equations numerically. For 
5 000 bins of the spectrum a system of 30 000 integro-differential equations describing the 
neutrino density evolution should be solved simultaneously. 

Hopefully, the asymmetry effect on CN for large mass differences may be successfully 
analytically estimated using the analytical approximations, accounting for the spectrum 



distortion of neutrino, like the one proposed in |45] for small mixing angles and large 
mass differences. Such approximations maybe very useful for analyzing oscillations effect 
on CN, because they simplify the equations and reduce the computational time. 

2.3.3 Relic lepton asymmetries. 

The case of different initial asymmetries, namely 10~ 10 < L < 10~ 4 , in models of CN 
with nonresonant oscillations, effective after electron neutrino freeze-out, was precisely 
studied in |St], [31], |3J], [J7| . It was found that depending on the concrete values of the relic 



asymmetry, it may suppress, enhance or not influence oscillations, and hence lead to an 
under-, overproduction of helium or to not changing its abundance at all. 

The possibility of lepton asymmetry to enhance oscillations, besides its well-known 
ability to suppress them, was revealed thanks to the study of neutrino evolution at each 
momentum: 



10 



For the parameters of the model the local term L dominates over the other terms for 
the mean neutrino momenta, and hence the resonant condition is not fulfilled: \L\ 3> \Q\, 
\L\ > cos 2$. 

However, for neutrinos with a given momentum p < p it is fulfilled and these neutrinos 
suffer a resonant transfer leading to the decrease of their number densities and a decrease of 
the asymmetry (since the resonance condition is fulfilled only in one sector, either neutrino 
or neutrino one, the particle densities in the other sector do not change). Then, due to 
the L decrease, resonant transfer becomes possible for neutrinos with higher momentum, 
and so on till L changes sign and the running resonance wave similarly passes through 
the antineutrino ensemble, leading to a rapid increase of L again till the next change of 
sign of L, etc. i.e. the observed enhancement is "spectrum resonance" effect. 

The total effect of the process is enhanced resonance transfer both in neutrino and 
antineutrino sectors, leading to greater overproduction of helium-4. Oscillation constraints 
for the case of CN with oscillations and with initial L were obtained [|3l[ . 



A similar investigation for the resonant case will be more complicated due to technical 
problems because it will deserve much greater number of bins for the spectrum distor- 
tion description because of the asymmetry growth, and hence, longer calculational time. 
However, this investigation is interesting, as far as such small values of the initial relic 
asymmetry are not excluded neither from observations nor from some profound theoretical 
principle. 

In conclusion we would like to stress that even a very small asymmetry, either initially 
present or dynamically generated, thanks to its indirect effect through neutrino oscillations 
is capable to influence strongly CN. Hence it should be accounted for as precisely as 
possible. 

2.4 Summary of neutrino oscillations effect on CN 

In case neutrino oscillations between active and sterile neutrinos proceed in the primordial 
plasma during the CN epoch, they can effect CN in the following ways: 

(a) In case active neutrinos have not decoupled, their oscillations to sterile ones may 
bring additional degrees of freedom into the primordial heat bath. This will lead to 
an increase of the Universe expansion rate H(t) and to an earlier freezing of the n/p - 
ratio, at times when neutrons were more abundant in comparison with the SBBN. Hence, 
this effect leads to an overproduction of helium-4. This effect was historically the first 
discussed |5| and on its basis first constraints on the oscillation parameters 5m 2 , sin 2 (2$) 
were estimated |g ^ ^ 0, 

(b) Oscillations considerably influence CN by distorting the neutrino and antineutrino 
spectrum. As far as the oscillation rate is energy dependent r osc ~ 5m 2 /E the low energy 
neutrinos start to oscillate first, and later the oscillations become noticeable for the more 
energetic neutrinos. Due to that, the spectrum of the neutrinos (antineutrinos) may 
become strongly distorted, especially in the case of oscillations into less abundant sterile 



neutrinos. This effect was shown considerable both in the vacuum oscillations case f|TU 



and matter oscillations case |29|. Spectrum distortion of the electron neutrinos due to 
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oscillations leads to overproduction of helium-4, which may be several times larger than 
He-4 overproduction due to an additional neutrino type. 

(c) Oscillations can lead to a substantial depletion of the number densities (N v and 
Np) of active neutrinos, in case they proceed between active and less abundant sterile 
neutrino states. This slows down the weak rates, T w ~ N V E%, and again leads to an 
earlier n/p - freezing and a corresponding increase of the helium-4 yield. 

This effect was analytically calculated and found to be important for CN first in the 
vacuum oscillation case JTU] . For matter oscillations with great mass differences it was first 
estimated in ref . |L7|, [2(J] , where it was described in terms of an effective chemical potential 
of neutrinos, generated due to active-sterile oscillations. Electron neutrino depletion was 
numerically calculated, without the account of spectrum distortion and asymmetry growth 



in ref. |[24| . Later these effects were taken into account in ref. |29[ for the case of small 
mass differences {8m 2 < 10~ 7 eV 2 ). 

For oscillations with small mass differences the latter two effects can lead up to 32% 



overproduction of helium-4 [BH . 

(d) Oscillations produce neutrino-antineutrino asymmetry, which on its turn influences 
the evolution of the neutrino and antineutrino ensembles and the oscillation pattern. In 
the nonresonant case the neutrino-mixing produced asymmetry was shown to have a 
negligible role on CN. However, in the resonant oscillation case the asymmetry effect 
on CN was shown to be considerable [[29], [33]]. For the case of small mass differences 
it was proven that even very small asymmetries L <C 0.01 considerably influence CN 
through oscillations, and therefore asymmetry effect on CN should be accounted for during 
asymmetry's full evolution. 

In general, dynamically produced asymmetry at small mixing angles suppresses os- 
cillations, which leads to less overproduction of helium-4 in comparison with the case 
of CN with oscillations but without the asymmetry account. Hence, the bounds on the 
oscillation parameters provided without the asymmetry account are alleviated at small 
mixing angles |3~3], |3~4 . 



Vice versa, the presence of a relic asymmetry in the nonresonant oscillation case leads 
to an alleviation of the bounds at large mixings due to suppression of oscillations by asym- 
metry. While at small mixings, due to spectrum resonance enhancement of oscillations 
caused by the asymmetry, the bounds on oscillation parameters are strengthened [31]. 

In case of neutrino oscillations, a precise description of the neutrino evolution at each 
momentum is necessary, inorder to account for the essential effects of spectrum distortion 
and asymmetry generation due to oscillations. To study the oscillation effects (a)-(d) 
on the light element production in CN, a self consistent study of neutrino and nucleons 
evolution is to be provided. 
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3 Updated constraints on oscillation parameters and 
the solar neutrino problem 



Inorder to account for all different effects of oscillations on CN we have made a selfcon- 
sistent numerical analysis of the kinetics of the oscillating neutrinos and the nucleons at 
freeze-out. 



3.1 The kinetics 

For simplicity we discussed a toy model of oscillations just in the electron sector z/j = Uu ui, 
I — e,s: 

v x = cv e + su s , 
v 2 = -sv e + cu s , 

where v s denotes the sterile electron antineutrino, c = cos($), s = sin(i?) and d is the 
mixing angle in the electron sector, the mass eigen-states v 1 and u 2 are Majorana particles 
with masses correspondingly m\ and m 2 . 

We have solved selfconsistently the set of the following coupled integro differential 
equations describing the evolution of the neutrino density matrix p and neutron number 
densities n n : 

Mt) _ R dpjt) | 
dt v dp„ 

+i [H , pit)} + iy/2G F (±£ - QIM 2 W ) iV 7 [a, pit)} , (1) 
(dn n /dt) = Hp n (dn n /dp n ) + 

+ J dQ(e~,p, v)\A(e~p -> un)\ 2 [n e -n p (l - p LL ) - n n p LL (l - n e -)\ 

-J dil(e + ,p, u)\A(e + n -> pu)\ 2 [n e +n n (l - p LL ) - n p p LL (l - n e +)\ . (2) 

where = U* e Uj e , p u is the momentum of electron neutrino, n stands for the number 
density of the interacting particles, dQ(i, j, k) is a phase space factor and A is the ampli- 
tude of the corresponding process. The sign plus in front of C corresponds to neutrino, 
while minus - to antineutrino. 

The initial condition for the neutrino ensembles in the interaction basis is assumed of 
the form: 

'1 0\ 
/ 



p = rC 



where n e u q = exp(-E v /T)/(l + exp(-E v /T)). 

7i is the free neutrino Hamiltonian. The 'nonlocal' term Q arises as an W/Z propa- 
gator effect, Q ~ E u T . £ is proportional to the fermion asymmetry of the plasma and 
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is essentially expressed through the neutrino asymmetries C ~ 2L Ue + L v + L Vt , where 
L^ T ~ (N^T—Nji^/Nj and L Ue ~ / d 3 p(pLL — Pll)/N 1 . Due to the different temperature 
dependence an interesting interplay between these two terms during the cooling of the 
Universe is observed. 

The neutron and proton number densities, used in the kinetic equations for neutrinos, 
are substituted from the numerical calculations of eq. (2). On the other hand, pn and Pll 
at each integration step of eq. (2) are taken from the simultaneously performed integration 
of the set of equations (1). 

These equations provide simultaneous account of the different competing processes, 
namely: neutrino oscillations, Hubble expansion and weak interaction processes. 

For the nonequilibrium active-sterile oscillations this is the only acceptable way of 
description of the asymmetry and oscillating neutrinos evolution. Because, in case of 
strongly distorted by oscillations neutrino and antineutrino spectrum, the asymmetry 
cannot be described in a degeneracy terms anymore, and the spectrum may strongly 
differ from an equilibrium spectrum with shifted effective temperature . 

The first equation results into a set of coupled nonlinear integro-differential equations 
with time dependent coefficients for the components of the density matrix of neutrinos. 
The number of these equations can be reduced to 6 equations for each momentum mode 
of neutrinos and antineutrinos, due to conservation of the total neutrino number density 
in the discussed model. The spectrum distortion was described by 1000 bins for the 
nonresonant case and by up to 5000 bins for the resonant case. In case the spectrum 
was described by N bins, a system of 6N + 1 coupled integro differential equations was 
numerically solved. 

The numerical analysis was provided for the characteristic temperature interval [2 MeV, 
0.3 MeV] [] and the full set of oscillation parameters of the active-sterile oscillation model 
|29| ], namely 5m 2 < 10~ 7 eV 2 and sin 2, $ > 0.001. We calculated precisely the n/p - freez- 
ing, which is the essential for the production of helium-4, till temperature 0.3 MeV, and 
accounted adiabatically for the decays of neutrons till the start of nuclear reactions below 
0.1 MeV. 

Updated v e <-> v s constraints 

Spectrum distortion, neutrino depletion and neutrino-antineutrino asymmetry of the os- 
cillating neutrinos, as well as the selfconsistent account of the evolution of neutrinos 
and nucleons, is essential for estimating oscillations effect on CN. In the last years the 
constraints in the electron-sterile sector were updated, accounting precisely for these ef- 
fects [0-0,0. 

The combined constraints for the nonresonant and the resonant case of electron-sterile 
oscillation parameters are shown in Fig. 2 for different values of relative increase of helium- 
4, 5Y p = (Y osc — Y p )/Y p . On the left-hand side of Fig. 2 the results for the nonresonant 
case are presented. The selfconsistent account of neutrinos and nucleons evolution and the 

6 At higher temperatures T > 2 MeV the deviations from the standard CN are negligible for the 
discussed model of oscillations with small mass differences. 
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precise calculation of neutrino depletion and spectrum distortion, allowed to strengthen 
the previous constraints by almost an order of magnitude at large mixings. 




Figure 2: The iso-helium contours corresponding to 3%, 5% and 7% overproduction of primor- 
dial helium abundance. LOW sterile solar solution is given by the closed dashed curves. 

The analytical fit to the updated constraints, corresponding to 5Y P = 3% is: 

<W(sin 2 2$) 4 < 1.5 x 1(T 9 eV 2 at 5m 2 > 0, 

5m 2 < 8.2 x KT 10 eV 2 at 5m 2 < at large mixing. 

In Fig. 3 the plot corresponding to Y p = 0.24 is compared with previous constraints: 
the pioneer estimates of refs. |17|, pOfl , partially accounting for neutrino depletion; the 



numerical calculations [ 23] , accounting partially for the neutrino depletion but neglecting 
the spectrum distortion and the dynamical asymmetry; the results of a recent analytical 
study P6 |, estimating the spectrum distortion effect in the approximation of small mixing 



angles and large mass differences. 

In the resonant case, due to the asymmetry growth account, He-4 overproduction is 
not so strongly expressed, and hence the old bounds on the oscillation parameters |24[ 



provided neglecting the asymmetry, were alleviated at small mixing angles (r.h.s. of 
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Figure 3: Update of the cosmological constraints for the electron-sterile oscillations. For com- 
parison, the dashed curves show the results from previous studies neglecting the spectrum dis- 
tortion and the growth of asymmetry [17, 20, 24], as well as the recent analytical constraints^] 
in the nonresonant case. The precise analysis, accounting for the spectrum distortion and the 
asymmetry growth [3C, 33] allows an almost complete exclusion of the LOW solution of the solar 
neutrino puzzle. 



Fig. 3). However, the precise account of the spectrum distortion of the oscillating neutrino, 
and the exact kinetic approach to both the neutrino evolution and to the nucleons freeze- 
out, strengthen the cosmological constraints at large mixing by an order of magnitude. 

It will be appropriate to provide similar investigations for large mass differences in 
the electron-sterile case and also to update the available constraints for the <-> v s and 
v T <-> v s cases by taking into account the energy spectrum distortion and the asymmetry 
growth. We expect that the constraints on the oscillation parameters for the <-» v s 
cases will be more slightly influenced than in the electron-sterile case, since fi and r 
neutrinos do not directly participate in the nucleon kinetics. 
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CN constraints on solar neutrino solutions 



These constraints on active-sterile neutrino oscillations exclude the active-sterile LOW 
solution to the solar neutrino puzzle in addition to the LMA solution, excluded in the 
pioneer works. LOW electron-sterile solution was obtained from the analysis of the 1258 
days SuperKamiokande experimental data on neutrino electron scattering and zenith angle 
variations of the solar neutrino flux |4T| , [43H . It is shown in Figs. 2 and 3 by the closed 
solid line around maximal mixing and 5m 2 ~ 10 -8 eV 2 . 

Assumed SY P < 3%, electron-sterile LOW solution is almost completely excluded for 
5m 2 < and it is completely excluded for 5m 2 < case. It is interesting to note also 
that even in case of very high primordial helium-4 SY P /Y P = 7%, sterile LOW solution 
still remains partially excluded [[37]. 

According to the global analysis of the solar neutrino data from SuperKamiokande, 
SAGE, GALLEX + GNO and the Chlorine experiments, the LMA and LOW solutions 
were found not acceptable for oscillations into sterile neutrinos |]40|, il, 421. The recent 



global analyses confirmed that LMA and LOW sterile solutions are disfavoured |44 



The conclusions from the global analysis of experimental data in 2000 ^0, ^2], |41| and 
in 2001 |43, H] are in remarkable agreement with the cosmological constraints on LMA 



(dating from the early 1980's) and on LOW solutions (obtained in 1999). 
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